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We have revealed that the barrier-tunnehng process in 
non-integrable systems is strongly hnked to chaos in complex 
phase space by investigating a simple scattering map model. 
The semiclassical wavefunction reproduces complicated fea- 
tures of tunneling perfectly and it enables us to solve all the 
reasons why those features appear in spite of absence of chaos 
on the real plane. Multi-generation structure of manifolds, 
which is the manifestation of complex-domain homoclinic en- 
tanglement created by complexified classical dynamics, allows 
a symbolic coding and it is used as a guiding principle to ex- 
tract dominant complex trajectories from all the semiclassical 
candidates. 
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Tunneling phenomenon is peculiar to quantum me- 
chanics and no counterparts exist in classical mechanics. 
Features of tunneling are nevertheless strongly influenced 
by the underlying classical dynamics In particular, 

chaotic features appearing in tunneling have been paid 
attention to in connection with real-domain chaos 

A promising approach to see the connection of these 
two opposite concepts is to carry out the complex semi- 
classical analysis, which allows us to describe and inter- 
pret the tunneling phenomenon in terms of complex clas- 
sical trajectories Q]. It has been shown that the complex 
semiclassical theory can successfully be applied even in 
classically chaotic systems and the origin of characteris- 
tic structures of the wavefunction inherent in chaotic sys- 
tems is explained by the complex classical dynamics [||. 
A significant role of almost real-domain homoclinic tra- 
jectories in the energy barrier tunneling has been pointed 
out based on the trace formula approach Q . Recently, it 
is found that fringed pattern appears in the wavefunction 
of the two-dimensional barrier tunneling problem as a re- 
sult of interference between oscillatory Lagrangian man- 
ifold 1^ . They have shown a detailed scenario describing 
how such interference emerges in accordance with the di- 
vergent movement of singularities on the complex t-plane. 

In the present paper, we shall report the strong con- 
nection between the barrier-tunneling process in non- 
integrable systems and the chaos in complex phase space 
by analyzing a simple scattering map model. In partic- 
ular, it will be shown that even though the real-domain 
classical dynamics exhibits no chaos, i.e. null topologi- 
cal entropy, complex-domain chaos can make tunneling 
process complicated. Moreover our present analysis sug- 
gests that chaotic tunneling can be understood in a uni- 



fied manner from the viewpoint of complex-domain chaos 
not only in case of dynamical but also energy-barrier tun- 
neling processes. 

We first introduce a scattering map model which is 
described by the following Hamiltonian 

= T{p) + V{q)Y,5{t^n), (la) 

n 

T(p)=pV2, V{q)^keM~iq^) (fc,7>0). (lb) 

A set of classical equations of motion is given as 
{q]+i,P]+i) = iqj,+T'{pj),Pj - V'{qj+i)), where prime 
denotes a differentiation with respect to the correspond- 
ing argument. Note that the real- valued dynamics of 
our scattering map does not create chaos in contrast to 
the maps defined on the bounded phase space. This 
is because the system has only single periodic orbit, 
{q,p) = (0,0), and thus the topological entropy of the 
system is null. As shown in the inset of Fig. |l|, the stable 
and unstable manifolds of the fixed point (0, 0), yV'*(0, 0) 
and W"(0,0), oscillate without creating homoclinic or 
heteroclinic intersections. Any manifold initially put on 
the real plane is fully stretched but not folded perfectly 
so that it leaves away to infinity along W"(0, 0). 

An incident wave packet is put in the asymptotic re- 
gion. The initial kinetic energy is given far less than the 
potential barrier located around the origin. The form of 
initial wave packet is given by 
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where C is the normalization constant, a is the squeezing 
parameter, and qa,Pa are configuration and momentum 
of the center of mass, respectively. 

|(g|C/"|^)| for n = 10 is shown in Fig. where U 
denotes the unitary operator of 1-step quantum propa- 
gation. Although the mean energy of the wavepacket is 
less than the barrier height, we can observe various struc- 
tures such as crossover of amplitude, existence of plateau 
regions, erratic oscillation on them, cliffs and so on. Sim- 
ilar structures are also found in case of dynamical tun- 
neling in the presence of chaos Q, and we hereafter call 
such characteristic structures "plateau-cliff structure" as 
a representative of typical features which are completely 
absent in one-dimensional tunneling. However, an essen- 
tial difference of the present situation from the previous 
one is that the plateau-cliff structure in the scattering 
model emerges under the situation where the real-valued 
classical dynamics creates no chaos. 
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In order to carry out complex semiclassical analy- 
sis, we first prepare a pair of new canonical variables 
as g = (g - ipa2)/(y2a), P = {p - iqa'^) / {^a~^) , 
and some notations as Qq = Qigo^Pa), Po = P{<1QiPq), 
Qa = Q{qa,pa), Pa = PiQajPa)- The Ti-stcp quantum 
propagator is represented as the n-fold multiple integral, 
and the saddle point evaluation is applied to yield the 
classical equations. Semiclassical Van Vleck's formula 
for K-step wavefunction takes a form as 
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where the sum is over complex classical orbits satisfy- 
ing the boundary condition, Im g„ = 0, which is neces- 
sary since we here observe our wavefunction as a func- 
tion of q„. A denotes the normalization factor, and 
S, (p are the action and the Maslov index of each con- 
tributing orbit respectively. W is the generating func- 
tion which gives canonical transformations such that 

dW/dq,, = p„, dW/dPo = -Qo- 

In order to represent the complex orbits which can con- 
tribute to the semiclassical propagator (H) , we introduce 
a variable, AQq = Qq — Qa- Using AQo the contributing 
orbits in (||) are represented as a set of initial conditions 
Mn = {AQo I AQo e C, Imq^iQa+AQo, Pa) = 0|. 

A typical pattern of Al„-set is shown in Fig. 0(a). 
Each string, which we call the branch hereafter, repre- 
sents an individual classical orbit appearing the semiclas- 
sical sum (5) and a single string covers the whole range 
(— cx), +oo) of the final q„. 

Though the morphology of the set is quite elusive, 
when sufficiently blowing up any local area with lots of 
branches accumulated, one can find the self-similar struc- 
ture schematically demonstrated in Fig. ^b). At the 
center, a chain-shaped structure is developed in the hori- 
zontal direction. A sequence of chains with finer scale are 
arranged in both sides of the previous chain, and does the 
same around any of those smaller chains, and so on. Then 
it may be natural to assign a notion of the generation to 
each chain as shown in the figure. Higher generations 
sprout from lower generations almost vertically, and we 
call these self-similar configurations the multi-generation 
structure(MGS). 

We can discuss significant branches embedded in the 
hierarchical structure, which determine the semiclassical 
wavefunction (j^. The amplitude of an individual con- 
tribution in the semiclassical sum is almost governed by 
Im5, an imaginary action, which implies that the orbit 
whose imaginary part is the smallest survives as the final 
semiclassical contribution. Furthermore, the amount of 
ImS is roughly estimated by how deep each trajectory 
passes through the complex domain. This means that 
the history of each trajectory in complex domain is a 
primary factor concerning its weight in the semiclassical 
sum. It implies that the semiclassical wavefunction could 
not be reproduced untill one finds how those orbits with 
the smallest imaginary actions behave, since it is almost 



impossible to extract significant trajectories out of a huge 
number of branches without any guiding principle. 

To this end, we begin with the precise definition of 
MGS, which is given by considering the stable and un- 
stable manifolds in complex phase space, in particular an- 
alytical continuation of the stable manifold W{0, 0). In 
order to connect MGS with >V*(0, 0) in complex space, 
we first introduce the normalized coordinate on W{0, 0). 
Let $ be the conjugation map from C to yV"*(0,0), sat- 
isfying the relation !F{£,) = = A^^^, where A 
denotes the largest eigenvalue with respect to the saddle 
point (0,0) 0. 

The normalized coordinate ^ is used in Fig. ||(a) 
and (b) to represent the intersection between >V''(0,0) 
and I, the initial value plane consisting of the whole 
(q,p)'s corresponding to {AQo I AQo G C}. Comparing 
both figures, one can recognize the intersection pattern 
has self-similarity. The intersection between >V"(0,0) 
and W{0, 0) shows almost the same self-similar pattern, 
which implies the homoclinic entanglement of the stable 
manifold in complex phase space This means that 
null topological entropy in the real domain does not ex- 
clude existence of chaos in the complex domain. 

The generation can be assigned to each point on X n 
yV*(0, 0). Suppose r be the minimum distance on ^-plane 
from (0, 0) to intersection points ^"^(J n >V^(0, 0)), and 
let D = G C I 1^1 < r}. Here the normalized coordi- 
nate plane is decomposed into disjointed annulus as, 

r''iD\:FiD))nr\D\T{D)) = (m^n) (4) 

\JP\D\HD)) = C (5) 



Then each annulus T''^{D\T{D)) plays a role specifying 
the individual generation. More precisely, if a point of 
I n W(0,0) is contained in $ • T~''[D\T[D)), we say 
the point belongs to the n-th generation. 

Next we show each point of X n 'W{Q, 0) can be asso- 
ciated with a single chain structure on A^„-set which is 
shown as the hatched zone in Fig. §(c). The final strings 
of those displayed in the figure after 10 iterations of our 
map are projected in real phase space in Fig. ^(d). The 
hatched part in Fig. ^(c) develops to the bold curves in 
Fig. |^(d), which almost coincide with the real-domain 
unstable manifold, yV"(0, 0) n 

This behavior is not limited to a particular chain- 
shaped structure, but every chain structure in A^„-set 
becomes close to W"(0, 0) n R^. Chain-shaped structure 
is always created around each point of X n W*(0, 0), the 
reason of which can be explained by describing the pro- 
cess of time evolution of a tiny domain containing a point 
of Xn W*(0, 0). Such a domain is guided to the real phase 
space by yV"'(0,0) and then smeared over yV"(0,0). De- 
tails will be reported in [Q. In this way, we can find 
one-to-one relation between a single chain structure in 
A4„-set and an intersection point of X n >V''(0, 0). Since, 
as explained above, the latter forms the self-similar or 
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generation structure reflecting the homoclinic entangle- 
ment of the stable manifold, MGS of Al„-set is controlled 
by that of the stable manifold in complex phase space. 

Since we have defined the generation in such a way, 
the orbits belonging to higher generations usually gain 
larger Im S, since, taking roundabout ways in the com- 
plex chaotic region, they approach the origin later than 
those of lower generations. In order to extract significant 
trajectories out of them, it is essential to use symbolic 
dynamics in complex phase space. 

In the present case, a nice symbolic space can be or- 
ganized as the union of a single-element set {O} and a 
direct product of the sign of Re q, Im q and a multiplicity 
index v E N. The symbol 'O' describes monotonic con- 
vergence to the origin (0, 0), and the index v appears as a 
reflection of the transcendental property of the potential 
function such that our initial-value problem has infinite 
number of solutions. Since frequent flipping of the sign 
or large modulus of g-component which is represented as 
large result in large imaginary actions, our final prin- 
ciple turns out to extract the sequences of symbols out of 
MGS which represent those orbits with no sign flipping 
and minimal q-components with 1^=1. It considerably 
reduces our task in searching significant complex orbits 
only to linear dependence of the time step, otherwise ex- 
ponentially diverges. More detailed explanations to con- 
struct symbolic dynamics in the complex-domain chaos 
will be reported elsewhere |ll| ]. In particular, complex 
Henon map is analyzed to elucidate the relation between 
the MGS and Juha set Hi. 

Fig. ^ displays the behaviors of intersection points 
in MGS as a function of time, together with the cod- 
ing sequences. Fig. ||(a) shows a typical behavior such 
that both Re q and Im q oscillate in an erratic manner 
for some initial time steps and eventually approach the 
origin. Fig. |4|(b) shows such an orbit trapped by com- 
plex period-2 orbits as evidence of regular behavior cm- 
bedded in MGS. Such itinerating or oscillating complex 
orbits contribute in principle, however, their amplitudes 
are much smaller than those of orbits with Re q and Im q 
decreasing monotonically to zero as shown in Fig. ^c). 

The semiclassical wavefunction evaluated finally takes 
the form in Fig. pL The agreement with the quantal 
one shown in Fig. ^ is excellent. The origin of plateau- 
cliff structure is understood by resolving the superposed 
semiclassical function into contributions from individ- 
ual generations. The plateau-cliff structure is created 
because each semiclassical component itself has a flat 
plateau accompanied by sharp drop and interference be- 
tween branches forming the generation give rise to erratic 
oscillation on the plateau. Discontinuity of amplitude, 
as found in g « —168, —207, —224, etc. are caused by 
the Stokes phenomenon which is inevitable in the saddle 
point method, and is treated in an appropriate way [p^ . 

It should be stressed that the most crucial step in our 
semiclassical analysis is to decode embedded information 
in MGS which reflects complex homoclinic entanglement. 
This means that emergence of the entanglement in com- 



plex phase space is an essential ingredient in our descrip- 
tion. We think it is a quite general event in chaotic sys- 
tems, and thus so-called plateau-cliff structure, which is 
typically observed pattern of tunneling in the presence of 
chaos 1^, must appear as a manifestation of such com- 
plex structures. We can therefore expect that, whether 
energetic or dynamical, there is a common semiclassical 
mechanism of the tunneling phenomena in chaotic map 
systems, which should be attributed to the structure of 
entanglement, typically observed as MGS. 
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FIG. 1. |(g|;7"|*)| for n = 10 calculated quan- 
tum-mechanically {h = l, cr^lO, fc = 500, 7 = 0.005, qc,~ —123, 
Pa =23). An incident wave packet is set in the side of g < 0. 
The center of mass has been already reflected by the poten- 
tial barrier at this time step. (Inset) W(0, 0) and W(0,0) 
merely oscillating without homoclinic entanglement. 
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FIG. 2. (a) A^„-set for n = 10. Those branches which we 
finally take into account for the reproduction of the wavefunc- 
tion are represented by the line styles such as dash, bold dash, 
bold solid, bold dot, bold dash dot and bold dash dot dot. 
Each line style corresponds to the one in Fig. ^ (b) Schematic 
representation for MGS. Blown-up figures of the second and 
third generations display the self similarity of MGS. Solid 
squares represent the intersection points with VV''(0, 0). 



(c) 




Re 5 




FIG. 3. (a) The intersection points belonging to the 
1st, 2nd or 3rd generation plotted on W(0,0). (b) Mag- 
nification of (a), displaying the 1st and 2nd generations, (c) 
Chain-shaped manifolds in A^„-set drawn schematically, (d) 
The final strings of those displayed in (c) after 10 iterations 
of our map projected in real phase space. The hatched part 
in (c) develops to the bold curves, which almost coincide with 
the real-domain unstable manifold, W"(0, 0) n R^. 
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FIG. 4. Behaviors of intersection points in MGS. Solid 
lines as Re q, and broken lines as Im q. The sequences of 
symbols coding them also attached. These are the trajectories 
(a) showing stochastic motions, (b) temporarily attracted by 
a complex period-2 orbit, (c) approaching real phase space 
monotonicaUy. 
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FIG. 5. The semiclassical wavefunction for 10 time steps 
resolved into contributions from different generations. Each 
component corresponds to the branches with the same line 
styles in Fig. ^(a). The thin dashed lines represent those 
branches which are dominated, in each coordinate q, by some 
other branches. (Inset) The semiclassical wavefunction finally 
obtained by the superposition of these contributions. 
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